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In this short review, the calculation of the next-to-next-to-leading order QCD corrections 
to the inclusive radiative decay B — > X„7 is described. I summarize the salient features 
of the calculational framework adopted, discuss the results obtained in the last few years, 
and indicate the technical tools that made the NNLO calculations possible. I conclude 
by comparing the current NNLO theoretical estimate for the branching ratio with the 
experimental measurement and by briefly discussing the size and origin of the residual 
theoretical uncertainty. 
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1. Framework 

The inclusive radiative decays of B mesons into a photon and an arbitrary hadronic 
state of total strangeness —l,B—> X s j, currently play a relevant role within the 
precision tests of the Standard Model (SM) and of its extensions. The decay process 
is sketched in Fig. [U B denotes a -B or B° mesons, while X s indicates an inclusive 
hadronic state not containing charmed particles. 

At the parton level, the process in Fig. [1] is induced by a Flavor Changing 
Neutral Current (FCNC) decay of the b quark contained in the B meson. The b 
quark decays into a strange quark plus other partons, collectively indicated by the 
symbol Xj arton , and a photon. In the SM, such a decay takes place at first at one 
loop, through "penguin" diagrams such as the one shown in Fig.O FCNC decays are 
rare (i. e. loop-suppressed) in the SM; therefore, they are very sensitive to Beyond 
the SM physics effects, which can arise in the perturbative expansion at the same 
order as the leading SM contribution. 

In contrast with the exclusive decay modes, inclusive decays of B mesons are 
theoretically clean observables; in fact, it is possible to prove that the decay width 
F(B — ► X s j) is well approximated by the partonic decay rate T(b — > X^ alton j): 

T(B -> X sl ) = T(b -> AT Ito " 7 ) + A non - port - . (1) 

The second term on the r. h. s. of Eq. (|TJ) represents non-pcrturbative corrections. 
The latter are small, since they are suppressed at least by a factor (AQCD/mb) 2 , 
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Fig. 1. A schematic illustration of the B — > Xs~f decay. 



where m& is the 6-quark mass and Aqcd ~ 200 MeV. The relation in Eq. fl} is 
known as Heavy Quark Expansion (for a review, see Ref. [TJ . 

The partonic process can be studied within the context of perturbative QCD. 
However, the first-order QCD corrections to the partonic process are very large. 
The large corrections originate from hard gluon exchanges between quark lines of 
the one-loop electroweak graphs (see Fig. [3]). In general, Feynman diagrams involv- 
ing different mass scales (say m\ and m,2, for example), depend on logarithms of 
the ratio of the masses. If there is a strong hierarchy among the two mass scales 
(e.g. mi -c rn-z), then the logarithms are numerically large. In the case of QCD 
corrections to the partonic process b — > Jf to °7, the mass scales involved are the 
VF-boson mass M w , the top quark mass mj, and mj,. M w and m t are of the same 
order of magnitude /i w ^100 GeV, while the 6-quark mass is considerably smaller: 
mi, ~ 5 GeV. At n-th order in the perturbative expansion in the strong coupling 
constant a s (mb), one finds terms of the form 



where m < n. For n = 1 and m = 1, the product in Eq. ^ is too large to be used 
as an expansion parameter; terms enhanced by large logarithmic coefficients must 
be resummed at all orders. Conventionally, calculations performed by resumming 
logarithms in which m — n arc referred to as leading order (LO) precision calcula- 
tions. By resumming terms in which m = n, n— 1, it is possible to obtain results 
of next-to-leading order (NLO) precision. Similarly, resumming m = n, n — 1, n — 2 
logarithms, it is possible to achieve next-to-next-to-leading order (NNLO) precision. 

The easiest way to implement the resummation of the large logarithms discussed 
above is to work within the context of a renormalization-group-improved effective 
theory with five active quarks. In such a theory, the heavy degrees of freedom 
involved in the decay under study are integrated out. By means of an operator 
product expansion, it is possible to factorize the contribution of the short-distance 
and long-distance dynamics in the decay of the B meson. In the SM, the short- 
distance dynamic is characterized by mass scales of the order of the top-quark or 
VF-boson mass, while the long-distance dynamic is characterized by the 6-quark 
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Fig. 2. One of the penguin diagrams contributing to b — ► sj decay. 

mass. The boundary between short-distance and long-distance is chosen at a low- 
energy scale fib such that rrib ~ fib "C M w . Clearly, the scale /Xf, is unphysical, and 
therefore physical quantities should not depend on it. However, all calculations are 
performed up to some fixed order in perturbation theory. This truncation of the 
perturbative series induces a dependence on the low-energy scale in the physical 
observables, which is formally of higher order with respect to the precision goal 
of the calculation. From the practical point of view, the Lagrangian employed in 
calculating the b — > X^' 10 "^ decay rate can be written a^f] 

C = £ QED®QCD ( u > d ,c,s,b) + 5^-^V t ;V r tt Ci(/i,/* w )Oi(/*) +® (jj^J • ( 3 ) 

In the equation above, £qjt;]3 and £qcd represent the usual QED and QCD 
Lagrangians with five active quark flavors. The second term in Eq. ([3]) is more 
interesting: Gf is the Fermi constant, Vt s and Va> are elements of the CKM matrix, 
and Oi are eight effective operators of dimensions five and six. Operators with 
dimensions larger than six are suppressed by inverse powers of the Vy-boson mass 
and are ignored. Finally, the short-distance dynamic is encoded in the "coupling 
constants" of the effective operators, which are called Wilson coefficients and are 
indicated by Ci in Eq. ([3]). The Wilson coefficients are the only elements of the 
Lagrangian which depend on the heavy particles masses M w and m t . 

Any perturbative calculation of the b — > Xf" 01 ^ decay rate within the context of 
the renormalization-group-improved perturbation theory applied to the Lagrangian 
in Eq. (J3j) requires three different steps: 

i) The first step, conventionally called matching, consists in fixing the value 
of the Wilson coefficients at the high-energy scale /i w ~ M w ,m t . This is 
achieved by requiring that Green functions calculated in the full SM and 
in the effective theory provide the same result up to terms suppressed by 
the ratio between the external momenta and M w or m t . At the scale /i w , 
QCD corrections are small and subsequently can be calculated in fixed-order 
perturbation theory 

ii) Secondly, once the value of the Wilson coefficient at the electroweak scale 
has been obtained from the matching step, it is then necessary to obtain 

a Terms suppressed by the CKM factor V* s V u i, are ignored; however, their NLO effect was ac- 
counted for in the NNLO calculation of the B — ► X s -y branching ratio in Ref. [2] 
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Fig. 3. One of the diagrams contributing to the first order QCD corrections to b — » s-y decay. 

the value of the Wilson coefficients at the low-energy scale fib ~ mj. This 
can be achieved by solving the system of renormalization group equations 
(RGE) satisfied by the Wilson coefficient, which have the following form: 

MX^G") =7j<0*)Cj(p)- ( 4 ) 

The matrix 7 in the equation above is the anomalous dimension matrix 
(ADM) of the effective operators. Since the various operators mix under 
renormalization, this step of the calculation is called mixing. By solving the 
RGE, it is possible to resum at all orders the large logarithms of the ratio 
jU w //ift in the Wilson coefficients. This becomes evident by looking at the 
leading order solution of Eq. ((4]) (for the diagonalized ADM, indicated by 
a tilde) 

-C/2/30 _ 

Ci(fi w ), (5) 

where Pq = 11 — 2/3 Nf is the leading order QCD beta function, and Nf is 
the number of active quarks, 
iii) Finally, it is necessary to calculate on-shell matrix elements of the partonic 
process in the effective theory. QCD radiative corrections to the matrix 
elements do not include large logarithms, since the dependence on the heavy 
degrees of freedom is completely encoded within the Wilson coefficients. 

The eight effective operators appearing in the Lagrangian in Eq. ([3]) are listed in 
|Appcndix A| At LO, the b — > amplitude is proportional to the effective Wilson 
coefficient C7 11 ; the effective coefficients were introduced in Ref. [31 

Radiative decays of the B meson were first experimentally observed (in the 
exclusive B — > K*^ decay mode) by the CLEO collaboration at Cornell in 1993. 
Nowadays, the branching ratio of the inclusive decay B — > X s j has been measured 
by several collaborations. The current world average, obtained by averaging the 
CLEO, BELLE, and BABAR measurements (Refs. H El U), is (Ref. [9]) 

B (B -> X s 7)^ >So = ( 3 - 52 ± 0-23 ± 0.09) x 10" 4 . (6) 

In Eq. ^ the first error is due to statistical and systematic uncertainty, while the 
second is due to theory input on the &-quark Fermi motion. In order to eliminate 
irreducible backgrounds, experimental collaborations impose a lower cut on the 
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photon energy. The value in Eq. ([6]) refers to a lower cut Eq = 1.6 GeV. In view of 
current experimental accuracy, theoretical predictions based upon NLO calculations 
in the a s expansion are no longer sufficient. In the last eight years, partial results 
necessary for the theoretical calculation of the B — > X s ^i branching ratio with NNLO 
precision were obtained by several groups (Refs.[I3E]EHE4E4E4[i6lElElEl 
[20J EJ)- Two years ago, by employing these partial results, it was possible to obtain 
the first theoretical estimate of the branching ratio at NNLO (f2ll2Tj). Currently, the 
NNLO program is heading toward completion (Ref . H21 [23l [24]) . 

In order to calculate the B — ► X s j branching ratio at NNLO, it is necessary to 
make use of some of the most powerful and recent techniques available for multi- 
loop Feynman diagrams calculations. The aim of the rest of this short review is 
to indicate the techniques employed in the parts of the NNLO calculation carried 
out so far, as well as to indicate the areas in which work is currently in progress. I 
conclude by comparing the current theoretical prediction of the branching ratio with 
the experimental measurements and by briefly discussing the residual theoretical 
uncertainties. 

2. NNLO Calculation 

The experimental error on the B — ■> X s j branching ratio in Eq. (|6|) is of 7%. 
The theoretical error affecting the NLO prediction is about 10% (see for example 
the comprehensive reviews in Refs. [23 [26]) . Aside for the NLO QCD corrections 
(Refs. [23E3C11E3E3E2H NLO determinations of 

the branching ratio include electroweak effects (Refs.l40" l l4 H l42 |) ; moreover, the nu- 
merical impact of the charm- and bottom-quark masses was analyzed in detail in 
Ref. H3j However, the NNLO QCD corrections were estimated to be at the level of 
±7% already in Ref. HU The inclusion of the NNLO QCD corrections is expected 
to significantly reduce the uncertainty associated to the charm-quark mass renor- 
malization scale. In fact, since the charm-quark mass m c first enters the branching 
ratio at NLO, the related scale dependence is a NNLO issue (sec Ref. 03]) • It is thus 
mandatory to include the numerically leading NNLO QCD effects in the theoretical 
prediction of the B — ► X s j branching ratio in order to bring the theoretical un- 
certainty at the same level of the experimental one. Below matching, mixing, and 
matrix element calculation at NNLO in QCD are discussed. 

2.1. Matching 

The goal of the matching procedure at NNLO is to calculate the order a 2 s corrections 
to the Wilson coefficients evaluated at the high-energy scale /i w ~ m t ,M w . The 
matching of the four-quark operators 0±, - - ,Oq, which requires the calculation 
of two-loop Feynman diagrams in the SM, was carried out in Ref. [TU) Since the 
magnetic- and chromo-magnetic dipole operators O7 and 0$ first arise at one-loop 
in the SM, the matching of these operators requires the calculation of three-loop 
Feynman diagrams, one of which is shown in Fig. 3] The dipole operator matching 




was carried out in Ref. [TTJ Both Ref. [TO] and Ref. [TTJ employ the same calculational 
technique. The SM diagrams are expanded in the ratio (external momenta) 2 //^ 
up to the first non-vanishing order. Spurious infrared divergencies originating from 
the expansion are regulated in dimensional regularization. All the masses, with 
the exception of M w , m t , and to&, are set equal to zero from the start; terms 
proportional to m 2 are also neglected. The integrals that must be evaluated after the 
expansion in the external momenta correspond to vacuum diagrams, involving one 
or, as in Ref. 111! two mass scales. Two- and three-loop vacuum integrals depending 
on a single mass scale are known. Three-loop vacuum diagrams depending on two 
different mass scales could not be evaluated directly. The problem was solved by 
expanding the integrals around the point m t = M w and for m t 3> M w . Both 
expansions give satisfactory results for the physical value M w ~ TOt/2. In order 
to write down the matching equations, it is necessary to require equivalency of 
the SM and effective theory Green functions. The calculation of the latter is not 
problematic; in fact, after expanding in the external momenta, loop diagrams give 
rise to scaleless integrals, which vanish in dimensional regularization. 



2.2. Mixing 

The entries of the ADM 7 in Eq. ([4]) can be obtained from the QCD renormalization 
constants in the effective theory. The latter are derived from ultraviolet divergencies 
in the Feynman diagrams with effective operator insertions. The ADM has the 
following pcrturbative expansion 

where the matrices at each order in a s have a block structure: 

/ A (k) R (k) \ 
7 (fe) = 7 X6 ■ (8) 

\ U 2 x6 ^2x2 / 

The lower-left block vanishes because the dimension-five dipole operators do not 
generate UV divergencies in dimension-six, four-quark amplitudes. The matrices 
and can be obtained from the UV divergencies of (k + l)-loop diagrams, 
while to evaluate it is necessary to evaluate the UV divergencies of (fc + 2)-loop 
diagrams. The matrices A^, B^ 2 \ and were calculated in the Refs. [T2J [TU and 
IT51 respectively. Sample diagrams needed in the calculation of the NNLO ADM arc 
shown in Fig. [5J Since one is interested in Feynman diagrams UV divergencies, it 
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(a) (b) (c) 

Fig. 5. Examples of Feynman diagrams needed for the calculation of the matrices (a), £?( 2 ' 

(b), and C< 2 > (c). Thick arrow lines indicate b quarks, thin arrow lines indicate massless quarks; 
the dark dots indicate effective operators. 

is necessary to evaluate the corresponding integrals by regulating UV and infrared 
singularities in different ways. In Ref s . [T2l [TH and [TBI this was done by introducing a 
common mass M for all fields, and thus by expanding the loop integrals in inverse 
powers of M. Therefore, the only integrals needed were single scale tadpoles up to 
four loops; these integrals are known. The number of diagrams involved is impres- 
sive: the calculation of alone requires the evaluation of over 20000 four-loop 
Feynman graphs. 

2.3. Matrix Elements - Total Rate 

To complete the calculation of the NNLO corrections to the b — ► X^ aTtan j decay 
rate, it is necessary to calculate the 0{a 2 s ) corrections to the matrix elements in the 
low-energy effective theory of Eq. ([3|). The decay width for the proces^] b — > X^ Ttoa "f 
can be written as 

r(b xr°"7K> Eo = G>c ^ )ml \v tb v;fY, °f^) cf 0*) G v M) • 

(9) 

The term proportional to Gij in Eq. ^ originates from the interference of diagrams 
mediated by the effective operator Oi and diagrams involving the effective operator 
Oj. Consequently, Gij is referred to as the (Oi,Oj) component of the decay width. 
In Eq. (JHJ) TOfc and mj, indicate the pole and running MS bottom quark mass, re- 
spectively. a em is the fine structure constant at zero momentum transfer. The total 
decay rate can be obtained from Eq. © by setting the lower cut on the photon 
energy, Eq, equal to zero. In this section we discuss the calculation of NNLO QCD 
corrections to the total decay rate. 

The (07,07) component is the only one completely known at NNLO in QCD. 
This component is numerically dominant and was independently calculated by two 
different groups (Refs. fTClfTTj) . In order to calculate the NNLO QCD corrections to 
the (07,07) component of the decay width in Eq. ([9]), it is necessary to consider 
three different sets of matrix elements: 

• Two-loop corrections to the process b — ► 37 interfered with the tree-level 

b I consider squared amplitudes summed over spin, color, and polarization of the final state, as well 
as averaged over the spin and color of the incoming b quark. 
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Fig. 6. Two-, three-, and four-particle cuts in a three-loop 6-quark self-energy diagram; they all 
contribute to the (O7, O7) component of the NNLO b — > X^' ton g decay rate. Thick lines indicate 
b quarks, thin lines indicate strange (massless) quarks. 



matrix element of the 6 — > S7 decay, 

• one-loop corrections to the process 6 — > sjg interfered with the correspond- 
ing tree-level amplitude, and 

• tree-level matrix elements for the processes b — > sjgg and b — > sjqq inter- 
fered among themselves. 

The contribution of each of the quantities listed above to the b — > Xj? arton 7 decay 
width is directly related to the imaginary part of three-loop 6-quark self-energy 
diagrams by means of the optical theorem (an example is shown in Fig. [H]). The 
contribution of each physical cut to the imaginary part of a self energy can be eval- 
uated by means of the Cutkosky rules f Refs. [TBI H6THT[) . In Ref. [T71 the contribution 
of each cut of each three-loop 6-quark self-energy diagram was evaluated by means 
of the methods commonly employed in the calculation of multi-loop Feynman di- 
agrams. By employing the Laporta algorithm (Refs. [48l |49l [50l [5"Tj) , it is possible 
to rewrite each cut diagram as a linear combination of a small number of Mas- 
ter Integrals (Mis). The Cutkosky rules associate a Dirac delta function with each 
line going through the cut. The latter sets the momentum flowing through a "cut 
propagator" on its mass-shell. The delta functions can be written as a difference of 
propagators (Refs. [52J, [53]); if q is the momentum flowing through a cut propagator 
one finds 

6 ( q2 + ma ) = - h ( g » + h + .0 ~ q- + h »o ) ■ (10) 

The reduction procedure is simplified by the fact that all of the integrals in which 
one of the cut propagators vanishes or is raised to a negative power are equal to zero, 
since in those cases the ±i0 prescription in Eq. (flO|) becomes irrelevant. In Ref. fTTl 
the Mis were evaluated numerically by means of the sector decomposition method 
(Ref. l54|) : the four-particle phase space integrals were parameterized according to 
the methods presented in Refs. l55l and 1561 After summing over the contribution of 
all of the cuts, the residual poles in the dimensional regulator e are canceled by 
UV renormalization. The methods employed in Ref. [IT] can also be applied to the 
calculation of the (Or, Os) and (0$, Og) components of the decay width in Eq. 
The authors of Ref. [TBI were able to analytically calculate the (6*7,07) component 
of the b — ► X^ Ita "'y decay width. In fact, in Ref. [57[ the same authors evaluated 
the complete imaginary part of individual three-loop 6-quark self energy diagrams. 
Each imaginary part is the sum of all of the possible cuts present in a given graph. 
In the case of the (Or,0$) and (Os,Os) decay width components, not all of the 
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cuts correspond to the interference of b — ► '"7 matrix elements. Therefore, in 
order to calculate the (07,0$) and (Os,Os) components, it will be necessary to 
individually evaluate each cut contributing to the process. 

The first calculation of b — ► x^' toI1 j matrix elements at NNLO was the one 
described in Ref. [15] In that work, a set of NNLO virtual corrections to the matrix 
elements involving the effective operator^ 0\,Oi, Or, Os, as well as bremsstrahlung 
corrections to matrix elements including Or, are evaluated. The authors restrict 
their calculation to diagrams with a closed light-quark loop. In fact, once the cor- 
rection proportional to oqNi (where Ni is the number of light quarks) are known, 
it is possible to estimate the complete NNLO corrections by means of the naive 
non-abelianization hypothesis (NNA, see Refs. l58l l59j) : the BLM (or large- /3q) ap- 
proximation is derived by replacing the factor Ni with the coefficient — 3/3n/2 in 
diagrams with a light-quark loop. In Ref. [HI the calculation was carried out by 
means of Feynman parameterization of the integrands, followed by an integration 
technique based upon Merlin-Barnes (MB) transform (see Refs. 1501 Ell EH)- The 
results are presented in analytic form. The three-loop virtual corrections to the 0\ 
and O2 matrix elements are expanded in powers of m^/m^. 

In Ref. [H Misiak and Steinhauser computed the NNLO corrections in the BLM 
approximation in the non-physical limit m c 3> m,b/2. They observed that their 
result, evaluated at m c ~ mt/2, matches well with the calculation of the same set 
of corrections in the small m c expansion (Ref. [T5]) . Because of this match, Misiak 
and Steinhauser also evaluated the complete m c -dependent NNLO corrections to the 
b — > x^ rton j matrix elements in the m c 3> m&/2 approximation. Subsequently, they 
assumed that the BLM result is a good approximation of the full NNLO corrections 
for m c — 0, and they interpolated their results for the non-BLM corrections down 
to the measured value of m c . The calculational technique employed in Ref. [2] is the 
same technique employed in the three-loop Wilson coefficient calculation presented 
in Ref. [TTJ The results of Refs. [15] and [2] were crucial to obtain the first NNLO 
estimate of the B — > X s ^f branching ratio (Refs 12] [21]). However, the calculation 
of the NNLO QCD corrections to the b — > Jf ,to °7 matrix elements is not yet 
complete, and several groups are still working to improve the understanding of this 
set of corrections. 

In Ref. EH the authors reconsider the virtual fermion loop corrections to the ma- 
trix elements involving the operators 0\ and Oi. They check the result of Ref. [TBI bv 
evaluating the corrections proportional to a^Ni, and they also numerically calculate 
the diagrams involving a closed massive charm or bottom quark loop. In the case 
of diagrams involving massless quark loops the reduction is carried out by means 
of the Laporta algorithm, while the Mis are evaluated in two different ways; first 
in an expansion in powers of m c /mb, and then by the numerical evaluation of the 

c At NNLO it is possible to neglect the operators O3, ■ ■ ■ , Oq, since they are suppressed by small 
Wilson coefficients; the NLO contribution to the branching ratio arising from these operators is 
< 1%. 
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Fig. 7. A sketch of the photon energy spectrum in the B — > X s -f decay. Only photons with energy 
Ey > Eo are observed. 

MB representation of the integrals. The Laporta algorithm is also employed for di- 
agrams involving a massive quark loop. However, for some of the Mis encountered 
in the latter case, the numerical evaluation of MB representations is not sufficiently 
precise. In these cases, an interesting technique based upon the numerical solution 
of the system of differential equations satisfied by the Mis is adopted. 

The effects of the charm-quark mass in the (6*7,07) and (Or,Os) decay width 
components are analyzed in Ref. [22] and Ref. [23] 

Currently, a group is working on the numerical calculation of the (0\,Oj) and 
(02,07) components in the m c = approximation (Refs. 163 ] [64 )) . The latter repre- 
sents the first step toward the calculation of the complete set of NNLO corrections 
in the m c = limit. Once this result becomes available, it will then be possible 
to interpolate the results of Ref. [2] to the physical value of m c without making 
any assumption on the behavior of the corrections for vanishing m c . In turn, this 
procedure will reduce the error associated with the m c interpolation. 

2.4. Matrix Elements - Spectrum 

In the two-body b — > 37 decay, the energy of the outgoing photon is fixed by 
the kinematic of the process: £7 7 = mt/2 in the 6-quark rest frame. However, two 
different phenomena give raise to a photon energy spectrum in the B — > X s j decay. 
First, the partonic decay in which one is interested is not b — > sj but rather b — > 
^parton-y. in events involving the emission of one or more gluon or qq pair, the photon 
has an energy Ej < m^jl. Secondly, the Fermi motion of the b quark in the B meson 
also contributes in generating a non-trivial photon energy spectrum. The situation is 
sketched in Fig [7] the smearing of the photon energy spectrum beyond the partonic 
end point mf,/2 is due to the non-perturbative Fermi motion effects, while the long 
low-energy tail has its origin in the gluon/quark-pair bremsstrahlung. The latter 
can be studied in pcrturbative QCD, while Fermi motion effects are modeled by 
means of a process- independent shape function (see Ref. [HI [55]) . Since experimental 
collaborations apply a lower cut on the photon energy of about 1.8 — 2 GeV, a 
detailed knowledge of the photon energy spectrum is mandatory in order to obtain a 
prediction for the B —> X s j branching ratio. The measured photon energy spectrum 
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provides direct information on the shape function, while from the moments of the 
truncated spectrum it is possible to extract relevant information on the Heavy 
Quark Expansion parameters. 

Working at NLO in a s , it is necessary to consider the virtual corrections to the 
process b — > 57 together with the matrix elements of the corresponding process with 
a gluon in the final state: b — » sjg. At NNLO, it is also necessary to consider the 
matrix elements with up to two gluons or a qq pair in the final state. After defining 
the dimensionless variable z = 2Ey/mb (0 < z < 1), the quantity Gij in Eq. ([9]) can 
be written as 



where one sets zq = 2Eo/mt, and /1 = mj,, and dGij/dz is the photon energy 
spectrum. 

The only component of the photon energy spectrum completely known at NNLO 
in QCD is G77 (Refs. [19l [20]) . which originates from the interference of graphs 
involving the magnetic dipole operator O7. This component of the spectrum is the 
numerically dominant one. The general structure of a spectrum component (Oi, Oj) 
is the following: 



Because of kinematic, only matrix elements with three or four particles in the final 
state contribute to the function R(z). Therefore, if the total decay rate Gij(zo = 0) 
is known, it is sufficient to calculate the function R{z) in Eq. (fT2|) . keeping z ^ 1, 
in order to know the photon energy spectrum. The constant / in Eq. (|12[) can be 
fixed a posteriori by requiring the equality of the total rate with the integral in 
Eq. (jlip once one sets Zq = 0. The (O7, O7) component of the total decay rate can 
be found in Refs . [TBI and 1 1 71 Consequently, in order to calculate the corresponding 
spectrum component, it is sufficient to consider one- loop corrections to the process 
b — ► sjg, as well as tree-level matrix elements contributing to the processes b — > sjgg 
and b — ► sjqq. At NLO, there are just two Feynman diagrams contributing to the 
(07,07) spectrum for z / 1. The squared amplitudes of the two graphs must be 
integrated over the final-state phase space, and the integrand must be multiplied 
by a Dirac delta function fixing the energy of the photon. Schematically, 



The optical theorem relates the decay width to the imaginary part of two-loop b- 
quarks self-energy diagrams. The contribution of each specific physical cut to the 
imaginary part of the two-loop sclf-cncrgy diagrams can be calculated by means of 
the Cutkosky rules. In particular, the quantity in Eq. (|13[) is obtained by summing 




(11) 



dGij (z) 



f ij 5(l-z) + R ij (z). 



(12) 



dz 
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Fig. 8. Cuts contributing to the b — > sg, b —* gg, and b — > sss processes in a three-loop b-quark 
self-energy diagram. 



over three different cut diagrams 
dG 77 (z) f /T\ 1 



dz NLO 

(14) 

where it is understood that the Dirac delta in Eq. (|13p was inserted by hand in the 
integrands. The same principle can be applied at NNLO, so that the (Oj, O7) photon 
energy spectrum can be obtained by calculating appropriate cuts of three-loop b- 
quark self-energy diagrams; examples are shown in Fig. [51 The contribution of each 
cut to the photon energy spectrum can be evaluated with the techniques discussed 
in the previous section. The Dirac delta function employed in order to fix the photon 
energy can be rewritten as a difference of propagators by means of the identity in 
Eq. (fT0|) already employed for the delta functions associated with cut propagators. 
The Mis can then be calculated in different ways. The differential equation method 
(Refs. [661 |67j) was employed in both Ref. Ql and Ref. US In Ref. H3 all of the Mis 
were checked numerically by means of the sector decomposition method. The sum of 
all of the cuts contributing to the magnetic dipole component of the spectrum is free 
from collinear and soft divergencies. In the intermediate stages of the calculation, 
the latter are regulated by dimensional regularization. The residual single pole in 
the dimensional regulator e is canceled by ultraviolet renormalization. The (O7, O7) 
spectrum is singular in the z — > 1 limit; the singularities appear in the spectrum 
as plus distributions of the form [ln"(l — z)/(l — z)] + (n = 0, . . . , 3). The singular 
terms were also predicted on the basis of the universality of soft and collinear gluon 
radiation (Ref. l68|) . It is interesting to observe that in the z — ► limit the (O7, O7) 
component of the spectrum should vanish like z 3 . This can be easily proved; each 
effective operator Oj gives raise to a factor z in the spectrum, and an additional 
factor z originates from the integration over the photon phase space. In a similar way, 
it is possible to determine the behavior of other components of the spectrum in the 
z — ► limit. Since the QED coupling of the quarks with the photon is proportional 
to 1/Ery for soft photons, the (Ot,Os) component of the photon energy spectrum 
is proportional to z when z — > 0. In the same limit, the (Os,Og,) component of 
the spectrum behaves like 1/z. These kinds of considerations are based upon the 
Feynman rules employed in the calculation; therefore, they apply to each cut of each 
diagram contributing to a given spectrum component. The qualitative information 
on the spectrum behavior in the z — > limit can be used as a powerful check of the 
calculation, but it can also be employed to fix some of the initial conditions needed 
to calculate the Mis with the differential equation method. In Refs. \W\ and [T71 the 
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c-quark mass is set to zero. The exact dependence of the NNLO corrections to the 
(O7, O7) spectrum and total decay rate on the c-quark mass is found in Ref. [221 

The techniques outlined above can in principle be used to also calculate the 
NNLO QCD corrections to the (Or,Os) component of the spectrum. In particular, 
the dependence of the (O7, 0$) spectrum on the c-quark mass was studied in Ref. 1231 

The BLM corrections to all of the components of the photon energy spectrum, 
with the exception of (Os,Og), were calculated almost a decade ago (see Ref. [T8| . 
It is well known that the BLM corrections, which are proportional to a 2 f3o, usu- 
ally provide a reliable estimate of the full a 2 corrections. In the case of the (O7, O7) 
spectrum component it was possible to verify that the BLM approximation provides 
a good estimate of the complete NNLO corrections. The BLM corrections to the 
(Os,Os) spectrum component were not calculated in Ref. IT8| they are suppressed 
by a factor (Qd.Cs/C7) 2 ~ 0.03, and therefore they are expected to be numerically 
small. Moreover, the (0$,Og,) spectrum component includes logarithmic singular- 
ities in the limit in which the strange quark is considered masslcss. As discussed 
above, the (Os, Os) spectrum is also singular in the limit of vanishing photon energy; 
therefore, its contribution to the branching ratio becomes numerically large for low 
(and experimentally unattainable) values of the cut on the photon energy. For com- 
pleteness, the BLM corrections to the (Os, Os) spectrum were recently computed in 
Ref. 1691 the collinear singularities were regulated by keeping a finite strange quark 
mass m s when needed. 

Finally, Table [T] summarizes the available NNLO QCD corrections to the matrix 
elements of the b — > X^'^j decay width, indicating both calculations of the total 
decay rate and of the photon energy spectrum. 

Tabic 1. NNLO QCD corrections to the various component of the b — > X J arton ^, decay width. 
The first row of the table refers to the full ct 2 . corrections, the second refers to the a, corrections 
with a light-quark loop, the third to a 2 corrections with a massive charm- or bottom-quark 
loop, and the fourth to non-BLM corrections in the m c 3> m^jl approximation. 



App. 


(Oi,Oi) 
(Oi,0 2 ) 

(o 2 ,o 2 ) 


(Oi,0 7 ) 

(o 2 ,o 7 ) 


(Oi,O a ) 
(0 2 ,Og) 


(0 7 ,0 7 ) 


(o 7 ,o 8 ) 
(o 8 ,o 8 ) 


full a 2 








Refs.ri6llT7llT9ll20l 




a 2 N t 


Ref.rrei 


Refs.[T8iri5ll24l 




Refs. [18] IT5l 


Refs. [T8iri5l [691 


c- b- loops 




Ref. [241 




Ref. EH 


Refs .1251 1551 


m c > rn b /2 


Ref. El 


Ref. El 


Ref. El 


Ref. El 





3. NNLO Estimate of the Branching Ratio 

According to the procedure of Ref. 1431 which is designed in order to reduce para- 
metric uncertainties originating from CKM angles and c- and 6-quark masses, the 
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B — > X s ^f branching ratio can be written as 



BIB 



X c ev) 



exp 



v t * s v tb 



VrJ 



6a P _ 



ttC 



[P(E ) + N(Eq)] , (15) 



where the perturbative corrections P(Eq) are defined as follows 



v cl 



6a e 



-P(E ) 



Vcl 



E>Eg 



T(b X u eP) 



and N(Eq) denotes non-perturbative corrections. The factor C is given by 



C 



cb 



r (g — > X c eD) 
r (B -> X„eP) 



(16) 



(17) 



Since C is a ratio of inclusive decay widths, it can be calculated by means of the 
same tools employed in the calculation of the B — > X s 7 width, and expressed as a 
double series in powers of a s and of Aqcd (see Ref . 1701 for a recent analysis). 
The current theoretical estimate of the quantity in Eq. (fT5|) in the SM is (Refs.[2I 



B (5 -> X s7 ) 



SM 



(18) 



^>i.eo.v= (3.15±0.23)xl0- 4 . 

The estimate in Eq. (fl8|) includes all of the numerically leading NNLO QCD cor- 
rections. When the value of C obtained in Ref. [70] is employed, the central value of 
Eq. (|18|) increases by a few percent. The error on the theoretical estimate is about 
7%, and was obtained by combining in quadrature four different uncertainties: para- 
metric uncertainty (3%), uncertainty due to missing higher order corrections (3%), 
uncertainty due to non-perturbative corrections (5%), and uncertainty due to the 
m c -interpolation ambiguity of Ref. [2] (3%). The estimate of Eq. (fT8]) does not in- 
clude some NNLO and non-perturbative corrections which are currently known; 
however, the combined effect of the neglected contributions is about +1.6%, which 
is a small correction compared to the theoretical uncertainty (more details can be 
found in Ref. [7T|) . It is important to observe that the result of Rcfs. [2] I2T1 reduced 
the theoretical uncertainty, which is approximately of the same magnitude as the 
experimental one. As expected, the inclusion of NNLO corrections significantly re- 
duced the dependence of the branching ratio on the matching scale fi w ~ A/ w , on 
the low-energy scale fib ~ mbi an d especially on the charm-mass MS renormaliza- 
tion scale /x c that first enters the calculation at NLO. The scale dependence of the 
LO, NLO, and NNLO predictions for the branching ratio are compared in Fig. 2 
of Ref [21] The SM theoretical prediction of Eq. ([18]) is now slightly more than la- 
lower than the experimental average in Eq. ([5]). Such an agreement can be used in 
order to set stringent constraints on the parameters of some Beyond the SM physics 
models (see for example Ref. |2"T]) . 

While progress in the calculation of perturbative corrections to the b — > x^ aTtan '-f 
decay width is expected in the future, the current theoretical error is dominated 
by the uncertainty associated to non-perturbative effects, estimated to be about 
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5% (Ref. [21]) . The non-perturbative uncertainty primarily arises from corrections 
of order 0(a s Aqcd/ti6) which are very difficult to evaluate; they were analyzed in 
Ref. [H 

In Table [21 I summarize the size of NLO, NNLO, and non-perturbative correc- 
tions to the branching ratio. 

Table 2. Size of various perturbative and non-perturbative set of corrections 
to B(B — > X s "y)^i 1 6GcV ^ ne percentages refer to the size of a given set 
of corrections with respect to a LO branching ratio of ~ 3.4 X 10 -4 . The 
corrections of 0(« s Ao_ciD/ m &) are n °t known. 



Perturbative 




Non-perturbative 




NLO QCD 0(a s ) 

NLO EW O(a) 
NNLO QCD 0(a 2 s ) 


30% 
4% 
10% 


LO QCD + NLO m b 0(A^ CD /m<;) 
LO QCD + NLO m c 0(A^ CD /m 2 c ) 
NLO QCD + NLO m b 0{a s A Q cn/m b ) 


1% 

3% 
5% 



In conclusion, the program which aims to calculate the NNLO QCD corrections 
to the b — > X^' toI1 j decay is well under way. The results so far obtained already have 
a substantial impact on the theoretical prediction of the B — > X s j SM branching 
ratio. The calculation of the NNLO corrections poses numerous technical challenges 
related to the number of Feynman diagram involved in various steps of the calcu- 
lation, as well as to the evaluation of the needed integrals. Such challenges require 
the extensive application of the most current calculational techniques developed for 
the automated evaluation of multi-loop Feynman diagrams. 

Appendix A. Effective Operators 

For completeness, the eight effective operators relevant for the b — > JSQ ,arton 7 decay 
are listed below: 

Ox 
3 

o 5 

6 
7 

o s 



= {s^T a P L c) {cj"T a P L b) , 
= (s^P L c) {cy"P L b) , 

= (s^P L b)J2(Qri) , 

9 

= {s ltl T a P L b)Y,^rT a q) , 
q 

= (si^ulpPLb) Y, irf-f-fo) > 

9 

= (si,l,l P T a P L b)J2(^Yl P T a q) , 
= 77T^ m b {sa^P R b) , 

= #2 TOh (^ UTap R b ) Ojo, . (A.l) 
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In the above equations, e and F^ v (g s and G^ u ) represent the electromagnetic 
(strong) coupling constant and field strength, respectively. The eight color genera- 
tors are indicated by T a . The sums in the operators O3, . . . , Oe run over the light 
quarks. 
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